Racunanje kompleksne struje generatora i snage na ulazu voda
Zc = sqrt((Rp + 1i*2*pi*f*Lp) / (Gp + 1i*2*pi*f*Cp))
Zc = 3.9673e+02 - 3.4586e+01i
gamma = sqrt((Rp + 1i*2*pi*f*Lp) * (Gp + 1i*2*pi*f*Cp))
gamma = 1.0832e-07 + 1.0706e-06i
a12 = Zc*sinh(gamma*D)
a12 = 5.0642e+01 + 3.2541e+02i
a21 = (1/Zc)*sinh(gamma*D)
U0 = a11*U + a12*ID
U0 = 1.4916e+05 + 1.6556e+05i
I0 = a21*U + a22*ID
I0 = 2.5169e+02 + 4.9322e+02i
S0 = U0*conj(I0)
S0 = 1.1920e+08 - 3.1899e+07i
Odredjivanje zamenske Pi sheme
syms Zc gamma D Z1 Z2 a11 a12 a21 a22 U0 I0 UD ID
assume(Z1~=0 & Z2~=0 & a11~=0 & a12~=0)
jednacine = [U0 == Z2*(ID+UD/Z1) + UD,...
I0 == U0/Z1 + (U0-UD)/Z2]
jednacine =

resenje = solve(jednacine, [U0,I0])
resenje =
U0: (UD*Z1 + UD*Z2 + ID*Z1*Z2)/Z1
I0: (2*UD*Z1 + UD*Z2 + ID*Z1^2 + ID*Z1*Z2)/Z1^2
jed1 = a11 == simplify(subs(resenje.U0/UD, ID,0))
jed1 =

jed2 = a12 == subs(resenje.U0/ID, UD,0)
jed2 = 
resenjeZ = solve([jed1, jed2], [Z1,Z2])
resenjeZ =
Z1: a12/(a11 - 1)
Z2: a12
Z1 = subs(resenjeZ.Z1, [a11 a12], [cosh(gamma*D) Zc*sinh(gamma*D)])
Z1 =

Z2 = subs(resenjeZ.Z2, [a11 a12], [cosh(gamma*D) Zc*sinh(gamma*D)])
Z2 = 